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Abstract 

Presented here is a calculation of inelastic J/ij: production in polarized 
photon-hadron collisions under the framework of NRQCD factorization for- 
malism. We consider the photoproduction of J /tp in the energy range relevant 
to HERA. The Weizsacker- Williams approximation is adopted in the evalua- 
tion of the cross sections for ep collisions. We found that this process can give 
another independent test for the color-octet mechanism, and the different fea- 
tures for the two color-octet processes may provide further informations on the 
mechanism for inelastic J/ip photoproduction. And the discrepancy on the 
production asymmetry A between various sets of polarized gluon distribution 
functions is also found to be distinctive. 



PACS number(s): 12.40.Nn, 13.85.Ni, 14.40. Gx 
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I. INTRODUCTION 



The nucleoli spin structure has attracted a lot of attention in the last ten years since 
the discovery of the so-called "spin crisis" by the European Muon Collaboration (EMC) 
During these times, much has been done to understand the spin structure of proton 
both on the theoretical side and the experimental side (mostly from the polarized DIS 
data). However, the use of DIS data alone does not allow an accurate determination for the 
polarized parton densities, especially for the polarized gluon distribution, which is crucial 
important to understand the gluon contribution to the proton spin. And now, the upgrade 
of HERA collider with polarized electron beam and proton beam at y/s = SOOGeV in the 
near future has been considered for a long time 0. With this high energy polarized electron- 
proton collider, one can investigate the polarized structure functions in much broader range 
of X and Q^, especially, to very low x and high region. And also the exploration of 
polarized gluon distribution AG(a;,(5^) will become accessible. To determine AG{x,Q'^), 
various photon-gluon fusion processes can be used. The typical promising process for this 
purpose is the large— di-jet production 0, which is a gluon induced process even in the 
leading order (LO) of QCD calculations. However, for this process one must deal with the 
uncertainty arising from the resolved photon contributions. In this paper, we will consider 
the inelastic J/ ijj production in the polarized photon-hadron collisions. Because the inelastic 
J/ip photoproduction is also a gluon induced process, so this process can be viewed as a 
possible way to study the polarized gluon distribution in the proton. Furthermore, the 
resolved photon contribution to this process is negligible in the major kinematical region 
of the produced J/ip. So, by inducing some kinematical cut, the uncertainties due to the 
resolved photon contribution can be eliminated safely. 

On the other hand, studies of heavy quarkonium production in high energy energies 
is another hot topic in recent years. In the conventional picture, the heavy quarkonium 
production is described in the color-singlet model [^. In this model, it is assumed that 
the heavy quark pair must be produced in a color-singlet state in short distance with the 
same angular-momentum quantum number as the charmonium meson which is eventually 
observed. However, with the recent Tevatron data on high px J/ip production appearing, 
this color-singlet picture for heavy quarkonium production became questionable. The data 
is larger than the theoretical prediction of the color-singlet model by a factor of about 
30 ~ 50 0. This is called the J/ip {ip') surplus problem. On the theoretical side, the 
naive color-singlet model has been supplanted by the nonrelativistic QCD (NRQCD) fac- 
torization formalism |p, which allows the infrared safe calculation of inclusive charmonium 
production and decay rates. In this approach, the production process is factorized into short 
and long distance parts, while the latter is associated with the nonperturbative matrix ele- 
ments of four-fermion operators. So, for heavy quarkonium production, the quark- ant iquark 
pair does not need to be in the color-singlet state in the short distance production stage, 
which is at the scale of l/mq {mq is the heavy quark mass). At this stage, it allows the 
other color configuration than the singlet for the heavy quark pair, such as color-octet. The 
later situation for heavy quarkonium production is called the color-octet mechanism. In 
this production mechanism, heavy quark- ant iquark pair is produced at short distances in a 
color-octet state, and then hadronizes into a final state quarkonium (physical state) nonper- 
turbatively. Provided this color-octet mechanism, one might provide explain the Tevatron 
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data on the surplus of J/ip and ip' production [^,0. 

Even though the color-octet mechanism has achieved some successes in describing the 
production and decay of heavy quarkonia, more tests of this mechanism are still needed. 
Recently, the photoproduction data from HERA puts a question about the universality 
of the color-octet matrix elements [|1^,|11|], in which the fitted values of the matrix elements 
{Ol^^ {}Sq)) and {O^^^ i^Pj)) are one order of magnitude smaller than those determined 
from the Tevatron data [0. (More recently, possible solutions for this problem have been 



suggested in |T2| , p!3| .) In this situation it is certainly helpful to find other processes to 



test the color-octet mechanism in the heavy quarkonium production. So, the inelastic J/ip 
production in polarized photon-hadron collisions can give another independent test of the 
production mechanism, aside from the unpolarized J/ip. 

The rest of the paper is organized as follows. In Sec. II, we will give the polarized cross 
section formula for the inelastic J/ip production at the polarized electron-proton collider. 
Here, we adopt the Weizsacker- Williams approximation to calculate the electroproduction 
cross section with the photoproduction cross section. The numerical results are given in 
Sec. HI. We will display the polarized cross sections for both the polarized photon-proton 
collisions and the electron-proton collisions. All of our results are relevant to the energy range 
of the upgrade polarized HERA collider. We also study the J/il) production asymmetry and 
its dependence on the polarized gluon distributions. In Sec. IV, we give the conclusions. 



II. POLARIZED CROSS SECTION FORMULAS 



In photoproduction processes, the electron-proton cross section can be equivalent to the 
photoproduction cross section convoluted by the photon flux factor from the electron under 



the Weizsacker- Williams approximation |14 
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Here, y is the energy fraction of the electron carried by the photon. /^'^•' is the flux factor of 
photon (or the distribution function of photon in the electron). In the following calculations, 



we use the form from 15 
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where rrie is the electron mass, and Q^ax is the maximum value of for the photopro- 
duction processes. In the following numerical calculations, we will adopt Q'^ax = IGeV'^. 
Compared with the original Weizsacker- Williams function, the above equation contains a 
non- logarithmic term with a singular behavior at y — |T^. This term will give non- 
negligible contributions to the unpolarized cross section. 

Correspondingly, the polarized electroproduction cross section has the similar form. 
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ep 



dyA&\y)dAa 
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(3) 



where the polarized photon flux factor function A/^^^ [y] has been derived in Ref. |T6| 
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This function also contains a non-logarithmic term, which however is not singular at y ^ 0. 
The numerical calculations show that this term will contribute about 5% to the total cross 
section of the inelastic J/ip photoproduction at ^/s = 300GeV for Q^^x = IGeV'^. The 
polarized cross section Aa in the above equation is defined in terms of (j{Xi, A2) with definite 
helicities for incoming particles, 

dAa = ^[da{+, +) + da{-, -) - d(T{+, -) - da{-, +)]. (5) 

In our study, we also adopt the Weizsacker- Williams approximation described above to 
evaluate the inelastic photoproduction of J/ip in electron-proton collisions. That is, the 
polarized and unpolarized cross sections for J/ip production in electron-proton collisions 
have the following forms, 

dAa{ep ^e + J/ij + X) = T™' dyAf\'\y)dAa{jp J/^j + X), (6) 

da{ep -^e + J/^ + X)= dy&\y)da{jp ^ J/^j + X), (7) 



where ymin and ymax are the kinematical boundary for the photon energy fraction from the 
electron beam. In the above equation, we do not include the resolved photon contribution 
to the J/xp photoproduction at ep colliders. This is because the resolved contribution is 
important only in the lower z region {z < 0.2) [jl^, where z is the energy fraction of the 
produced J/V". So, by inducing a lower bound cut on z, the resolved photon contribution can 
be eliminated. Aaljp —>■ J/ip + X) is the polarized cross section for the inelastic production 
of J/ip in the polarized photon-proton fusion processes. 

At high energies, the photoproduction processes proceed predominantly through photon- 
gluon fusion. So, the polarized cross section Aa{'jp J/ijj + X) is related to the polarized 
gluon distributions in the proton, 

dAa{-fp ^ J/iIj + X) = J dxAG{x, Q^)dAa{-fg J/tP + X). (8) 
The polarized gluon distribution function AG{x, Q"^) is defined as usual way, and now there 



are various parametrizations for this quantity, such the GRVS sets |18], the GS sets |19], the 
DSS sets [^. 



In photon-gluon fusion processes, according to the NRQCD factorization formalism, the 
unpolarized cross section and the polarized cross section can be factorized as the following 
forms, 

da{-f + g^J/^ + X) = J2F{^ + g^n + X) {O'J^), (9) 

n 

dAa{^ + g^J/ij + X) = Y: AF(7 + g^n + X){0'J^). (10) 

n 

Here, n denotes the cc pair configuration in the intermediate state (including angular mo- 
mentum '^^'^^Lj and color index 1 or 8). F{'~f + g —>■ n + X) and AF{'~f + g —>■ n + X) are 
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the short distance coefficients for the unpolarized and polarized cases respectively for the 
subprocess 'j + g — >■ n + X. {O^^'^) is the long distance non-perturbative matrix element 
which represents the probability of the cc pair in n configuration evolving into the physical 
state J /if). The short distance coefficients F and AF can be calculated by using perturbative 
QCD in expansion of powers of Us- The long distance matrix elements are still not available 
from the ffist principles at present. However, the relative importance of the contributions 
from different terms can be estimated by using the NRQCD velocity scaling rules. 

In order to suppress the diffiactive contribution and higher-twist corrections, we induced 
the lower bound cut for the J/V' transverse momentum > IGeV. With this constraints 
(non-zero px for J/ip), there are the following 2^2 partonic processes contributing to J/ip 
inelastic photoproduction in the photon-gluon fusion processes, 

'y + g^ccCSuD+g, (11) 

^ + g-^ccCSo,8)+g, (12) 

-f + g-^cc('S,,8)+g, (13) 

j + g^cc('Pj,8)+g. (14) 

The ffist subprocess is the color-singlet process, and the last three processes are all the color- 
octet processes. The unpolarized cross sections for these processes have been calculated in 
literature P JT0|JT7[1 . The polarized cross section for the color-singlet process has also been 
calculated in pT| . 

To calculate the polarized cross sections for the color-octet processes, we employ the 
helicity amplitude method. With this method, one can calculate the cross sections (j{Xi, A2) 



with definite helicities for the incident particles. Following p2[, we choose the polarization 
vectors for the incident photon, the incident gluon and the outgoing gluon as the following 
forms, 



= N,y^ ^2 ^3(1 T 75)+ ^3 ^2 ± 75)], (15) 
= N,[^2 h ^1(1 T 75)+ i>x i>z ^2(1 ± 75)], (16) 

= iVe[^3 ^1 ^2(1 T 75)+ ^2 #1 ^3(1 ± 75)]. (17) 



Where pi, p2, Ps and ei, 62, 63 are the momenta and the polarization vectors for the incident 
photon, incident gluon and outgoing gluon respectively. And the normalization factor is 

Ne = ^i==, (18) 

where the Mandelstam invariants s, t, u are defined as 

s = (pl+P2)^ i={.P2~Pzf, u={pi-p'if, (19) 
and they satisfy the relation 



s + t + u = M=Am 



With these definitions for the polarization vectors of the photon and gluons Eqs. (|TBp- 
(p!7|), the calculations of the helicity amplitudes are straightforward. From these helicity 
amplitudes, we can easily get the cross sections with the definite helicities for the incident 
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particles. And then we can obtain the polarized cross sections dAa{''yg ^ J/ip + X) for 



the different partonic processes Eqs. ([T^) -(14) from the definition Eq. @. Their analytic 



expressions are listed in the Appendix. The unpolarized cross section for every partonic 
processes can also be obtained by summing up the cross sections with all of the possible 
helicity states for the external particles. We have checked that our results for the unpolarized 
cross sections are consistent with the results of Ref. |T0|JT7[] . For comparison, we also list the 
unpolarized cross sections in the Appendix. 

After getting the polarized cross section, we can also calculate the asymmetry 

where j3 represents for some kinematical parameters, such as J/ip transverse momentum 
Pt and the energy fraction z. z is defined a.s z = p ■ kj/^jp ■ pi with p, kj/^, pi being 
the momenta of the proton, the outgoing J/ip and the incident photon respectively. The 
production asymmetry parameter A is more efficiency than the polarized cross section Act 
for the study of the polarized processes in the polarized collisions. This is because A is 
normalized, and then less depends on the input parameters which will strongly affect the 
normalization of the cross sections, such as the charm quark mass, strong coupling, and the 
factorization and renormalization scales. 



III. NUMERICAL RESULTS 

Provided the polarized and the unpolarized cross sections for the subprocesses, the in- 
elastic J/ip production rate and asymmetry in polarized photon-hadron collisions can be 
obtained. For the numerical evaluation, we choose the strong coupling constant a^, the 
charm quark mass mc, and the factorization scale /i^ to be 

as = 0.3, rric = 1.5 GeV, /i^ = i2m^f. 

For the parton distribution function of the proton, we use the GRV LO parametrization 
| 23| . For the polarized gluon distribution function, we will consider the GRVS sets [jl8[, the 



GS sets [ID, and the DDS sets || . 



We first display the different partonic processes contributions to the polarized cross 
section of J/ip inelastic production, where we adopt the GRVS STD set as the default 
set for the polarized gluon distribution function of the proton. Fig. 1 is for the photon- 
proton collisions, where we typically set the energy of photon-proton system to be W.yp = 
lOOGeV relevant to the photoproduction processes at HERA. Fig. 2 is for the electron- 
proton collisions, where we set the energy range also relevant to HERA, ^/s = SOOGeV^. For 
the electroproduction cross sections, the integral region of W^p is set to be: {25GeV)'^ < 
W^p < (ISOGeV^)^ typically used in the photoproduction processes at HERA. In these two 
figures, we both plot the z and pf. distributions: (a) are z distributions and (b) are p"^ 
distributions. To evaluate these curves, we have imposed a cut on z and P^: 0.1 < z for the 

distributions and > IGeV"^ for the z distributions. For the NRQCD long distance 
matrix elements, we set as 
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(OfCSi)) = l.lGGeV^ (21) 

{OtCSi)) = 1.06 X 10-^GeV\ (22) 

{OtCSo)) = 3.0 X 10-^GeV\ (23) 

{OtCPo))/ml = 1.0 X 10-^GeV\ (24) 

{OtePj)) = {2J+l){OtePo)). (25) 

The last equation comes from the heavy quark spin symmetry of NRQCD. The color-singlet 
matrix element (Of^^Si)) can be related to the wave function at origin, and can be taken 
its value from the leptonic decay width of J/ip. The value of the color-octet matrix element 
{Og(^Si)) is taken from a fit to the large pt J/ip production at the Tevatron |^|. This 



matrix element is not important to J/^jj photoproduction both for the unpolarized and 
polarized cases. On the other hand, the two other color-octet matrix elements, ((9g (^Sq)) 
and (Og (^Po)); are known to be very important in the inelastic J /"(p photoproduction p!0| , p!7| . 
However, their values are not well determined from the present experimental data on J /if) 



productions. Here, we just follow Ref. and take their values tentatively as listed above 
(which are also consistent with the naive NRQCD velocity scaling rules) to see what their 
contributions to the polarized cross section of the inelastic J/%Ij photoproduction will be. 
The solid curves in these two figures are for the color-singlet process contributions, the 
dashed lines for the color-octet ^5*0 contributions, the dotted lines for the color-octet ^Pj 
contributions, and the dotted-dashed lines for the color-octet ^5*1 contributions. 

From Figs. 1 and 2, we can see that in the lower z region of the inelastic J/ip production 
the dominant contribution comes from the color-singlet process. In the larger z region, 
the dominant comes from the color-octet ^5*0 process, and its contribution rise rapidly with 
z in the whole region. This property is as the same as that for the unpolarized cross 
section |]TD|. However, for the unpolarized cross section of inelastic J/ip photoproduction. 



we know that this rapidly rising property at large z is not consistent with the experimental 



measurements fg]. It is argued in [|17| that in this region, the soft gluon resummation is 
important, which may smear the z distribution in the large z region. So, the predictions for 
the LO calculations may be not rehable. Correspondingly, the predictions of the polarized 
cross section for the color-octet ^5*0 contributions (the dashed lines) of Figs. 1 and 2 in the 
large z region will also be strongly affected and modified by the higher order corrections. 
For the color-octet ^Pj processes, on the other hand, their contributions are always smaller 
than the color-octet ^5*0 contributions. And especially, their contributions at large z are 
similar to the color-singlet contributions, i.e., fall down rapidly with z at the end point 
region. This property is quite different from their contributions to the unpolarized cross 
section of J /ip photoproduction |1To| , [T7| , which also rise rapidly at large z like the color-octet 



^5*0 contributions. The p\ distributions of Figs. 1 and 2 are always dominated by the color- 
octet ^5*0 process contributions. And their contributions are larger than the color-singlet 
contributions by almost an order of magnitude in the whole region of which is similar 
to their contributions to the unpolarized cross section [|l0|,|l7l. However, the color-octet 



^Pj contributions to the p\ distributions of the polarized cross sections again have different 
features compared to their contributions to the unpolarized cross sections ||T0|JT7[| . Their 
contributions are much smaller than the color-octet ^5*0 contributions in the whole region 
of py, and even smaller than the color-singlet contributions at low At very low p^ 
(p|. ~ IGeV'^), their contributions to the polarized cross section Act even give a negative 
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value, which is not plotted in the figures. 

In Figs. 3 and 4, we display the asymmetry of Eq. ( pO]) for J / ip inelastic photoproduction 
as functions of z and p"^. Fig. 3 is for the photon-proton collisions, and Fig. 4 is for the 
electron-proton collisions. The solid curves are the production asymmetries in the color- 
singlet model. The other two curves represent the production asymmetries after including 
the color-octet contributions in two cases of values for the color-octet matrix elements: (I) 
for the dashed lines, {OtCSo)) = 3.0 x IQ-^GeV^ and {Ot{^Po))/m^^ = 0; (II) for the dotted 
lines, {OtCPo))/ml = 1.0 x IQ-^ GeV^ and {Ot{^ So)) = 0. The three curves in each diagram 
have the same behaviors as each other. For the z distributions, the asymmetries for all of 
the curves rise with z, and then fall down after reaching their maximum values. There is no 
steep rising as in the polarized and unpolarized cross section distributions when z approaches 
its end point for all curves. In the whole region of z, the case (I) color-octet predictions 
are larger than the color-singlet predictions, and the case (II) color-octet predictions on the 
other hand are similar to the color-singlet predictions in the lower z region and even smaller 
in the larger z region. For the distributions, the asymmetries for the three curves all rise 
with p"^. The case (I) predictions are similar to the color-singlet predictions. However, the 
case (II) predictions are smaller than the color-singlet predictions by almost a factor of two. 
Especially, in the low p"^ {p\ ~ 1 — 2GeV'^) region, the production asymmetry A for case 
(II) is very small comparable with zero. 

Concluding the results displayed in the above four figures (Figs. 1-4), we can see that the 
difference between the color-singlet and color-octet predictions on the distributions behaviors 
and the absolute sizes of the polarized cross sections Act and the production asymmetries A 
may provide some information on J/ip production mechanism. That is to say, the inelastic 
J /if) photoproduction in polarized collisions will give another independent test for the color- 
octet mechanism. In particular, as shown in the above analysis, the different features of the 
contributions from the two different color-octet processes may provide further important 
informations on the mechanism for the inelastic photoproduction of J/^. 

Comparing the results of Figs. 3 and 4, we can see that the production asymmetries for 
the polarized electron-proton collisions are all smaller than those for the polarized photon- 
proton collisions by about one order of magnitude. This is mostly because the polarized fiux 
factor A/^^) in Eq. (D has a minus sign in the coefficient before the logarithmic term (the 
dominant term) compared with the unpolarized flux factor /^*^) in Eq. (^, which will result 
in a big difference between them when convoluted with the photoproduction cross section 
by using Eqs. (|^) and (|^. 

Now, we turn to study the dependence of our results on the different parameterizations 
for the polarized gluon distributions in the proton. In Fig. 5 and 6, we plot the production 
asymmetry of J / ip as functions of z and p^ with different sets of the polarized gluon distri- 
bution functions [ll8H20[] (for all these sets, we just use the LO results of them). We have 
tried all of the polarized gluon distribution parametrizations for every sets mentioned above, 
but found the different parametrizations for each set give the similar results. However, the 
difference between different sets are much larger. For representation, we choose the polarized 
gluon distribution function for each set as follows: the GRVS STD parametrization for the 
GRVS stes (the solid lines), GS-C for the GS sets (the dashed lines), and DSSl for the DSS 
sets (the dotted lines). Fig. 5 are the z distributions and Fig. 6 are the p"^ distributions: 
(a) are for the results only from the color-singlet contribution; (b) are for the results after 
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including the color-octet contributions, where we choose the color-octet matrix elements of 

{Oti'So)) and {Oti'Po)) as 

{OtCSo)) = {OtePo))/ml = 0.008. 

From these two figures, we can see that the difference of the asymmetry A between the 
different polarized gluon distribution functions is distinctive both for the color-singlet and 
color-octet predictions. The GRVS STD set give the production asymmetry much larger 
than those for the other two sets, and the GS-C set and the DDSl set give similar results 
for the production asymmetry. 

IV. CONCLUSIONS 

We have calculated in this paper the inelastic J/ ip production in polarized photon-hadron 
collisions. We have evaluated the polarized production rate and the production asymmetry 
at polarized electron-proton colliders relevant to the HERA energy range. We performed 
our results both for the color-singlet model predictions and the predictions after including 
the color-octet contributions. We found that the inelastic J/ ip photoproduction in polarized 
collisions may give another independent test for the color-octet mechanism. Especially, 
the different features of the contributions from the two different color-octet processes may 
provide further important informations on the mechanism for the inelastic photoproduction 
of J/ip. And most important, we found that the discrepancy on the production asymmetry 
A between the various sets of the polarized gluon distribution function is distinctive. So, 
the process studied in this paper may be a candidate process for the measurement of the 
polarized gluon distribution function in the proton at future polarized HERA collider. 
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APPENDIX: 



In this appendix, we list the polarized and unpolarized cross sections for the different 
partonic processes of Eq. (|lT])-(|^). For convenience, we define two variables, P = st+tu+sit 
and Q = {§ + i){s + u){i + u). 

7 + (7 — > cc(^Si, 1) + g: the unpolarized cross section, 

da _ 2MiATTfaalel P' - M^stu ,^^^^,^ 

the polarized cross section, 

dAa _ 2M{Anfaalelius^ - P 



dt 277rs2 g2 



(orr^i)) (A2) 



7 + 5^ ^ cci^Si, 8) + g: the unpolarized and polarized cross sections are relevant to those for 
the process 7 + (7 — ^ cci^Si, 1) + g multiplied by the factor 



8 (Of(35i)) 

1 + g ^ cc(^S'o,8) +g: 



-{Ot{'P,))[^^ + 



2 



dt AttPQ^ N «V ^ M3(t + s)2(t + M)2 

^ M3t(t + S)2(s + U)2 M3t(t + S)2(s + U)2 ^' 

-f + g^cc{^Pi,8)+g: 

da (4:n)^aa'ief.{Ot(^Pi)),^r,^. ^„ „ ^ „ , „ 

= 27rg2Q4M3 + + + + ^ ) + ^ + ^) 

(3P + 7Pu + 7tu^ - u^) + s^(3P + 22t^M + GOt^w^ + 76pM^ + 36?m^ 
+Aiu^ - u^) + + u){t^ + 12pM + + 72Pm3 + 32Pm4 + Atu^ 



(A3) 



rfor _ ^jA-nfaayju M« + + + ix" ^^^^^^ ^ 

rf^ - ma^ MQ^ ^ ^^^^ 

-IT- m^^ Mq^ {o,[s,)). (A5) 

7 + ^ cc(3Po, 8) + (7: 

da {AT^faalel , .^^ ,3 „ , , .OM^gu su1?{2M^ + 3tM2 + gM)^ 



dt AttPQ^ n »v ^ M3(t + g)2(t + M)2 

^ M3t(t + g)2(g + M)2 ^ M3t(t + g)2(g + M)2 J' ^ ^ 

c/Aa {Anfaalel .^^^ ,3 ^ , , .gM^gw gMp(2M4 + 3tM2 + su] 



(A7) 
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+u^) + 2sHu{i^ + IQt^u + 38t^M^ + 59Pu^ + 38Pu^ + lOiu^ + u^) 

(2P + 9£m + 2^2) + i%^{t + (A8) 

^— = '-^pA-^:^-^\2s\t^ + tu + u^) + 2s^t t + u)i2t + 3u 

+sH{2P + 5Pu + 2iu^ + u^) - 2s^u{i + m)(2? + 8Pu + Stu^ + ii^) 

-siu{i + + 7tu + u^) - Pu'^{i + uf] (A9) 

7 + (7 cc(3P2, 8) + g: 

^ V — I s c\ 8V i2s\(t + + 12s^u(t + m St^ + 7tu + Au^) + / 

+ 142t^M + 38At\^ + 454t\^ + 303tM^ + 96m^) + + M)(9t^ + 202t^M 
+438t^u^ + 442t\^ + 309tM^ + 120u^) + s\9f + 198t^u + 736t^u^ + UOOt^u^ 
+1184t=^M^ + 860t^M^ + 429tM'' + 96m^) + + u){3f + I00t\ + A50t\'^ 
+720t^u^ + 688t\^ + A96t\^ + 255tu^ + 48m^) + 2sM'^lt^ + H^fu + 3621%"^ 
+585t^u^ + eOOt^u^ + 440t^u^ + 227t\^ + 66tu^ + 6u^) + 2sV{t + uf{l9t' 
+76t^u + im^u^ + 8At^u^ + 48tu^ + I2u^) + st^u^{t + uf{22t'^ + 59^3^ 
+58t\2 ^ 3g^^3 ^ ^2m^) + 3(t + u)H^u% (AlO) 

= ^ — " " \ — —\l2s^{t + u)2 + 12s^(t + u){At^ + 5tu + 3u^) + 2s^(41t^ 

+99^3^ + mt^u^ + 78tu3 + 30m^) + 6s^(t + M)(12t^ + 2lt^u + 24t2u2 + 20tM^ 
+10m^) + s^(30t'' + 117t^M + 22Qt\^ + 285t^M^ + 228t^M^ + imu^ + 24m^) 
+2sh{t + M)2(2t^ + 16t\ + 20t\^ + 13tu^ + 6u^) + + M)^(7t^ + 

-17tM2 - 12^3) + 3tV(t + n)3]. (All) 

In practice, the three "^Pf^ {J = 0, 1, 2) subprocesses contributions to J/i/j production 
are always added together, by using the relation of their color-octet matrix elements due 
to the heavy quark spin symmetry Eq. (^). For convenience, we also give here the total 
sum of the cross sections for the three ^Pf^ processes, which is more simpler than the above 
expressions. 

The sum of the three ^Pf'^ subprocesses, 

dai^pf) 3(ATTfaayjOt(^P2)), 7 ^ « . -^9 - 9^ 

+s^{i + u){2P + A5Pu + mu^ + 35u^) + s^(4p + 63i% + 132Pu^ 
+ 156?M^ + 98tM^ + 35u^) + s^i + m)(2P + A5i% + 91?^^ + 99Pu^ 
+57iu^ + 21u^) + s'^u{l3i^ + 70Pm + I36i%'^ + I32pu^ + 88?^^ 
+38tM^ + 7u^) + siu\i + n)(13t^ + 34Pn + 29Pu^ + 18iu^ + 7u^) 
+2i%^{i + u)% (A12) 
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(llP - 17tM + Sm^) + s^(9t^ - 28Fm - 48?^^ - 26tM^ + S-u^) 
+s^{i + u){8i^ - 19Pu - 37Pu^ - 23iu^ + 2u'^) + sH{2i + u) 

-ISiu^ - 10#) + 2i%\i + uy]. (A13) 
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Figure Captions 



FIG. 1. The polarized cross section distribution as a function of (a) z and (b) pf. in polarized 
photon-proton coUisions at the energy range of W^^ — lOOGeF. 

FIG. 2. The polarized cross section distribution as a function of (a) z and (b) in polarized 
electron-proton collisions at the energy range relevant to HERA, = SOOGeV^, where we 
use the Weizsacker- Williams approximation to evaluate the production rate for electron- 
proton collisions. 

FIG. 3. The production asymmetry distributions for Fig. 1. 
FIG. 4. The production asymmetry distributions for Fig. 2. 

FIG. 5. The z distributions of the production asymmetry A for the inelastic J/-?/' photopro- 
duction in polarized electron- proton collisions for different sets of polarized gluon distribution 
in the proton, (a) are the results for the color-singlet predictions, and (b) are the results 
after including the color-octet contributions. 
FIG. 6. The same as Fig. 5, but for the distributions. 
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TP collisions, E^^=100GeV, GRVS STD 
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